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The possibility of multiferroicity arising from charge or- 
dering in LuFe204 and structurally related rare earth fer- 
rites is reviewed. Recent experimental work on macro- 
scopic indications of ferroelectricity and microscopic de- 
termination of coupled spin and charge order indicates 
that this scenario does not hold. Understanding the ori- 
gin of the experimentally observed charge and spin or- 
der will require further theoretical work. Other aspects 
of recent research in these materials, such as geometrical 
frustration effects, possible electric-field-induced transi- 
tions, or orbital order are also briefly treated. 




Bilayer subunit in i?Fe204 compounds with proposed chai'ge 
order rendering it polar. 
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1 Introduction The interplay of magnetic and elec- 
tric degrees of freedom in correlated electron materials is 
one of the major topics of contemporary condensed matter 
physics, with a rapidly rising number of publications per 
year fT). Materials combining magnetism with ferroelec- 
tricity, termed multiferroics, offer the prospect of a large 
magnetoelectric coupling, which has a high potential for 
applications in future information technology [2|. Because 
the traditional mechanism of ferroelectricity is incompati- 
ble with magnetism |3!|, much interest has focused on un- 
conventional routes to ferroelectricity (see H for a recent 
extensive review). Particularly intriguing is the mechanism 
of ferroelectricity originating from charge ordering (CO) 
(reviewed in |I5]|6|), the ordered arrangement of different 
valence states of an ion, typically a transition metal. Be- 
cause any ferroelectric polarization is built from electric 
dipole moments, i.e. non-centrosymmetric charge distribu- 
tions, any CO breaking inversion symmetry automatically 
induces a polarization, which may be very large [iSj. The 
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presence of different valence states of a transition metal ion 
implies an active spin degree of freedom on the same ion, 
and therefore a strong magnetoelectric coupling may be ex- 
pected as well. This mechanism of multiferroicity is thus 
very attractive from the point of view of prospective ap- 
plications. However, while the concept is straight-forward, 
examples of oxides where this mechanism is experimen- 
tally indicated to occur are exceedingly rare and none is 
really well understood. 



The often cited (e.g. lIHISB) prototypical example ma- 
terial providing "proof of principle" for this mechanism 
is LuFe204, a charge- (and spin-) frustrated system. In 
the following, I will review recent research on this and 
isostructural materials, which may serve also as an illumi- 
nating case study for (much less investigated) other poten- 
tial examples of CO-based ferroelectricity. LuFe204 had 
been proposed in 2005 to be a multiferroic due to Fe^+/'^+ 
CO ifTTl . based on both macroscopic indications by dielec- 
tric spectroscopy and pyroelectric current measurements, 
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Figure 1 a) RSm Crystal structure of 
-RFe204 and description with different cells 
(after (7]). Apart from the primitive rhom- 
bohedral and the _R-centered hexagonal 
cells the lattice could also be described by 
one of three C-centered monoclinic cells, 
which are rotated by 120° with respect to 
each other. Cells describing the charge and 
spin order in the three domains are obtained 
from these by tripling the monoclinic b and 
doubling the monoclinic c axes (omitted for 
clarity), b) R^^ ion size effect on cell vol- 
ume and ratio of intralayer Fe— Fe distance 
ahcx to bilayer thickness db (c.f. c), com- 
piled from |8||9I10II11I12II13I14I . c) Iron 
bilayers with one of the monoclinic cells, 
and the corresponding supercell. Blue tri- 
angle illustrates geometrical frustration ef- 
fects (see text). Red lines denote a minimal 
set of interactions I15II16I necessary for 3D 
charge order. 



and on a plausible microscopic model of a polar CO. This 
proposed "ferroelectricity from CO" is the main focus of 
the review, in particular recent findings by structure re- 
finement of a non-polar CO IfTSl and the suggestion of the 
absence of polar order by dielectric spectroscopy HI 911201 . 
Other aspects that have been of recent focus in these mate- 
rials, such as geometrical frustration expressed in CO and 
magnetic order, strong spin-charge coupling, or potential 
electric-field-induced phase transitions, are treated as well. 
The main experimental findings reviewed are summarized 
in Sec. |71 which concludes with a brief outlook on future 
research directions both in this family of compounds and 
for "ferroelectricity from CO" in general. 

2 The i?Fe204 family Rare earth ferrites i?Fe204 
with R a 3+ ion without partially occupied d levels ( Y, Ho, 
Er, Tm, Yb, Lu, or In), known since the 1970s ||2T1 . crys- 
tallize in a rhombohedral (i?3m) structure ll8]l9] [T0lfmfT2l 
[T3lfT4ll (see Fig. [T^) featuring characteristic triangular bi- 
layers of trigonal-bipyramidal-coordinated Fe as the elec- 
tronically active subunit (c.f. Fig.[T];). Magnetism of the Fe 
spins attracted an initial wave of research into these materi- 
als, with unusual features such as anomalous thermomag- 
netization |22]| or giant coercivity |231| discovered. In all 
investigated compounds, a very large Ising anisotropy with 
spins pointing perpendicular to the layers was observed 
12311241 . The average valence of the Fe ions is 2.5+, and 
therefore Fe^"'"/'^+ valence order or CO may also be ex- 
pected, first investigated by Mossbauer spectroscopy ||251 . 
Hence, there are two binary degrees of freedom at the Fe 
sites; Ising spin t / ^ and valence 2 + /3+. 

Geometrical Frustration The arrangement of Fe 
ions leads to "geometrical frustration" ||261 hampering 
both spin order (SO) and CO, illustrated in Fig. [T];) for 
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Figure 2 Magnetization 
A/(r) of polycrystalline 
YFe204-5 with differ- 
ent S, measured in 0.4 T 
upon heating, after the 
sample has been cooled 
in the same field (o) or in 
zero-iield (•). Two tran- 
sitions in the most sto- 
ichiometric sample are 
marked. Data from 11271 . 



the Ising spins by a blue triangle: repulsive or antiferro- 
magnetic (AF) nearest-neighbor (NN) interactions within 
the layers can be satisfied only for two of the spin pairs 
of the triangle. For an isolated layer, this would lead to 
a macroscopic degeneracy of ground states. The proxim- 
ity of another triangular layer changes the topology of 
the problem, particularly since the inter-layer Fe-Fe NN- 
distance (marked U2 in Fig. [T]:) is about 8% shorter than 
the intra-layer one (marked Ui). However, there is poten- 
tial for frustration in the inter-layer interaction as well, 
because each Fe ion has three NN in the other layer (see 
also Fig.[8]l. Furthermore, the rhombohedral stacking leads 
to frustration of the interactions between different bilayers, 
hindering full 3D order Even when macroscopic ground- 
state degeneracy is broken, geometrical frustration often 
Il26l leads to a competition of different phases with nearly 
the same energy, facile creation of defects, and complex 
unusual ground states, such as the proposed (VJ\ ferroelec- 
tric CO in these materials, discussed in Sec. |4] 

Strong impact of oxygen stoichiometry A soon 
observed experimental difficulty in these compounds was 
a strong variation of physical properties for different sam- 
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Figure 3 Magnetization M{T) of three different LuFe204 
crystals measured in 100 Oe||c, upon cooling except where 
noted. Panel a) after 1281 , panel b) after | 



pies. For YFe204_5 this was clearly linked to tiny vari- 
ations in oxygen-stoichiometry ll27l . as shown in Fig. |2] 
magnetization M{T) for the most off-stoichiometric sam- 
ple shows a large difference between zero-field-cooled and 
field-cooled M suggesting "glassy" magnetism without 
establishment of long-range SO, whereas AI{T) for the 
stoichiometric sample indicates two distinct well-ordered 
antiferromagnetic phases. Single-crystal neutron diffrac- 
tion confirmed a tendency towards 3D long-range order 
for samples closer to stoichiometry [30), but due to the ab- 
sence of sufficiently stoichiometric crystals until recently 
BTI . the magnetic structures remain unsolved. Such a 
sample-dependence was observed also for Er ll32l and Lu 
(see Fig. |3]l, although the link to O-stoichiometry ll33l is 
less clear, mostly due to the range of stoichiometries for 
which long-range order is observable shrinking with de- 
creasing _R'^+ size (c.f. Il32l and Fig.[TJ)). In a wider range 
of stoichiometries, links to the remanent magnetization 
and resistivity are clear, however ll34l . 

Ion-size effects and structural modifications 
Changing the i?^+ size leads to an expected change of the 
cell volume (^ 5% between In and Y), but more relevant 
for CO and SO is a 15% change in the ratio between the 
layer-separation within a bilayer, db, and the intra-layer 
Fe-Fe distance ahox (Fig. [TJi), which can be expected to 
modify the relative importance of interactions Ui and C/2. 
Partial substitution, e.g. by Co, Cu, or Mn, on the Fe site 
generally leads to decreased conductivity, but also sup- 
pressed magnetic ordering (see e.g. [35. 136]). However, 
effects due to substitution can be difficult to disentangle 
from effects due to a potentially changed O-stoichiometry. 
Finally, a larger structural modification can be done by 
intercalating i?Fe204 with one or more blocks of i?Fe03, 
each adding an additional R— and a single-Fe layer be- 
tween the bilayers ll37l . There have been relatively few 
studies on these intercalated compounds (e.g. Il38l ). Fi- 
nally, a high-pressure orthorhombic modification with a 
large supercell, possibly CO, has recently been reported 
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Figure 4 Real part of the dielectric permittivity e' vs 
T (lines) in LuFe204 for frequencies between 1 Hz and 
1 GHz equally spaced with two frequencies per decade. 
The data additionally marked with symbols (0,+) are both 
taken at 100 kHz, but with electrodes from different mate- 
rial: Graphite (o) and silver (+ and all other data). Figure 
from [19 1 , @ 2012 American Physical Society. 



but the details of this structure have yet to be eluci- 
dated. 

3 Anomalous dielectric dispersion and pyroelec- 
tric currents Many rare earth ferrites, e.g. ErFe204 BOl . 
have been studied extensively by dielectric spectroscopy, 
with anomalously large real parts of the dielectric constant 
e' in the low-frequency limit generally observed. A typi- 
cal (recent) example of such a measurement is shown in 
Fig- 111 From the frequency-dependence of the temperature 
of maximum rise of £'{T) a connection of the dispersion 
to electrons hopping between Fe^+ and Fe'^+ has been 
concluded and an origin of the large dielectric constants 
in motion of ferroelectric domain boundaries suggested 
111711401 . Within this interpretation, giant (> 20%) room- 
temperature magneto-dielectric response iHTI would then 
suggest multiferroicity with large magnetoelectric cou- 
pling, potentially useful for applications. Research on 
possible ferroelectricity in rare earth ferrites focused on 
LuFe204, because in this compound superstructure reflec- 
tions were found, which for example by resonant x-ray 
diffraction at the Fe /v — edge could be clearly associated 
with Fe^+Z'^^ charge order 01711421 . allowing the discus- 
sion of potentially polar CO phases (see Sec.|4|i. 

Macroscopic proof of ferroelectricity requires demon- 
stration of a remanent electric polarization, which should 
be switchable by application of an electric field. A rema- 
nent polarization larger than in the traditional ferroelectric 
BaTiOa was suggested by pyroelectric current measure- 
ments ITtI . shown in Fig.|5^). These measurements were 
performed by first cooling the single-crystal sample in an 
electric field of ±10kV/cm, then switching off the elec- 
tric field, and finally measuring, upon warming in £' = 0, 
the pyroelectric current associated with the change of the 
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Figure 5 a) Electric polarization P vs temperature T 
of LuFe204 as deduced from pyroelectric current mea- 
surements conducted on warming after the sample had 
been cooled in electric fields E\\c of ±10kV/cm. Panel 
from ini, © 2005 Nature Publishing Group, b) P vs 
E hysteresis loops of LuFe204 at 140 K for two differ- 
ent frequencies, from |20J, © EDP Sciences, Societa Ital- 
iana di Fisica, Springer- Verlag 2012. With kind permis- 
sion of The European Physical Journal (EPJ). c) P{E) of 
Lu2Fe2.14Mno.86O7 at 2 kHz, from ||38l, © 2009 Ameri- 
can Institute of Physics. 



polarization. The dependence of the pyroelectric current 
and thus the polarization on the direction of the cooling 
field is consistent with a macroscopic polarization, and the 
decay to zero around the charge ordering temperature of 
Tco ^ 320 - 330 K El,, 44] supports a connection of P 
with the charge ordering, although the lack of a saturation 
region with no pyrocurrent flowing anymore above Tco 
renders this last conclusion tentative. It was noted, how- 
ever, that similar behavior can also be observed in (non- 
ferroelectric) leaky dielectrics, where it originates from lo- 
calization of free carriers at interfaces ll45l . 

Subsequent measurements by various groups partly re- 
produced the pyrocurrent results of Fig. [5^), as noted e.g. 
in II42I . or tried to establish P{E) hysteresis loops at low 
temperatures to more directly show ferroelectricity. Out of 
many attempts at the latter, only few were published Il20l 
[46" 47 1, some showing hysteresis, but not the saturation 
necessary [48j to establish polarization switching. Fig.|5}') 
shows P{E) measured at 140 K and two frequencies Il20l : 
the low-frequency loop shows substantial hysteresis, but no 
indications of saturation, whereas the loop measured with 
a frequency of 80 kHz exhibits practically linear behavior 
without appreciable hysteresis. 

Because the polarization is linked with the dielectric 
constant, P = (e — l)6oE, detailed modeling of the di- 
electric dispersion can show the origin of observed polar- 
ization behaviors. This was done for the results obtained 
on LuFe204 single crystals over 9 decades of frequencies 
|fT9l and shown in Fig. |4] The qualitative behavior with 
large e' at low frequency and high T is the same as in pre- 
vious observations flTl, but a strong impact of the elec- 
trode material as shown for the example of 100 kHz (differ- 
ent symbols) already suggests significant contact effects. 



e.g. due to Schottky-type depletion layers at the contact in- 
terfaces, which act as an additional capacitance lfT9l . The 
whole data-set of e'{f, T) and e"(/, T) could be weU de- 
scribed with an equivalent-circuit model featuring contact- 
capacitance and -conductance, and sample hopping- and 
dc-conductivity, in addition to the intrinsic dielectric con- 
stant £i of the material, all parameters independent of fre- 
quency. The fitted £i{T) is only around 30 and does not 
show any indications of ferroelectric or antiferroelectric 
transitions. Two recent studies on polycrystalline LuFe204 
reached the same conclusion 12011491 . For low frequencies 
or high T, the circuit reduces to a leaky contact capacitance 
in series with the sample resistivity, resulting in large e' . 
Because the step in e' mainly depends on contact capaci- 
tance and sample resistance |[T9l a magneto-resistance will 
lead to an apparent magnetodielectric response, according 
to the analysis in P^ the latter can be fully accounted for 
by the former. A very recent study |l50l of dielectric spec- 
troscopy using various contact materials suggests a small 
(~ 1%) sample magnetodielectric effect at least at 220 K 
(at the AFM-fM transition, c.f. Fig. [Tl]i, but the simulta- 
neously observed 18 times larger change in resistivity may 
influence the fitted sample capacity. 

P{E) loops can only show the intrinsic behavior if 
measured at frequencies and temperatures where e'{f,T) 
is close to £i{T). For the measurements shown in Fig.|5}'), 
this is the case for P{E) measured at 80 kHz, but clearly 
not for P{E) measured at 0.2 Hz. There is one P{E) loop 
in the literature |!38| that does seem to indicate ferroelec- 
tricity given the expected tendency to saturate in high E 
(Fig. |5};), measured on intercalated Lu2Fe307, in which 
conductivity has been suppressed by ^ 29% Mn for Fe sub- 
stitution. The occurrence of at least piezoelectricity in this 
compound is supported by piezoresponse force microscopy 
lISTI . Although the switchable polarization is very weak 
(compare panels c and a), smaller even than in "spin-spiral 
multiferroics" like TbMnOs, this suggests that intercalated 
compounds deserve more attention than they currently get. 

4 Charge Order It was the combination of macro- 
scopic indications of ferroelectricity, discussed in Sec. [3] 
with a likely microscopic model of charge order (CO) that 
is electrically polar, which made a convincing case for 
LuFe204 exhibiting ferroelectricity from CO IITtI . This 
CO-model was deduced by scattering methods, in which 
signatures of the order in reciprocal space are investigated. 

For many i?Fe204, at least diffuse scattering of x-rays 
at the ikk^) line in reciprocal space (in hexagonal-cell no- 
tation) suggested a tendency towards a CO that could be 
described, for example, with a so-called -\/3 x v^ cell. Al- 
ternatively, a propagation vector with {\\) in-plane com- 
ponent can be conveniently described by expressing the 
crystal structure without CO in a C— centered monoclinic 
cell, which is then enlarged three times along its 6— axis. 
Because such a CO breaks the 3 -fold rotation symmetry 
of the crystal structure, domains of CO with symmetry- 
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Figure 6 Sketch of the 4 symmetry-allowed Fe"^+ (black) 
/ Fe''+ (white) charge configurations following (|-|0) or 
(■ii|) propagation (with additional contribution by (000) 
or (00|), respectively). Arrows indicate bilayer polariza- 
tion. 



equivalent (||) and (||) propagations and 120° rotated 
monoclinic cells, as sketched in Fig.[Tl also exist 1441 . 

The simplest arrangement of Fe'^+ and Fe'^^ ions in 
a bilayer leading to such propagations and maintaining 
charge neutrality is the one proposed in ifTTI and shown 
in Fig. |6] top-left (two cells of the CO are shown). Here, 

.3+ 



High-resolution x-ray diffraction Ill5ll43ll44l found the 
exact positions of the (i, -i ,halfinteger) reflections to be 
well described by s + pa, s + pe, and s + pc, where 
s is a structural reflection following the rhombohedral 
centering condition —h + k + i = 3n with n integer. 
The three symmetry-equivalent propagation vectors, cor- 
responding to CO domains, are pa = {^ + 6,^ + S,^), 

Pb = (I - 2(5, 1 + 6, 1), and pc = (| + (5, | - 2S, |) and 
S ^ 0.003 T-dependent gH. A secondary set of reflec- 
tions is described by p^ = (r, r, |), pg = (2r, r, |), and 
p^ = (t, 2t, |), and t ^ 96 following the T— dependence 
of 6. Although this seems to suggest an incommensu- 
rate modulation of the Fe valence state, there is no evi- 
dence of a "true" incommensuration, which would imply 
a wide distribution of Fe valence states from 2+ to 3+, 
in conflict with spectroscopic measurements Il25ll63ll64l . 
Much more likely is therefore a locally commensurate 
state, interspersed with discommensurations or anti-phase- 
boundaries (Fig. [8]l of an average separation, estimated 
from T ~ 0.028 ll44l . of about 12 mii. It has been demon- 
strated on other charge-ordering oxides that a proliferation 
of such anti-phase-boundaries, which can cost very little 
energy due to geometrical frustration, can lead to sharp su- 
perstructure reflections at incommensurate positions |l 65l . 
Consequently, the possible CO corresponding to ( ^ -H) 
and for completeness also C|jO) have to be examined 
through symmetry-analysis II66II. For each of these prop- 
agation vectors, there are two irreducible representations, 
corresponding to the two Fe ions of the primitive unit cell 
(labeled 1 and 2 in Fig.[T^) having the same (J2) or differ- 
ent (Fi) valence. Each corresponding CO contains more 
than than two valence states, in contradiction to the re- 
sults of Mossbauer spectroscopy, which imply a bimodal 
valence distribution ll63l . Such a bimodal distribution is 
obtained by combining, in a unique way ll44ll . the CO from 



The four resulting CO are shown in Fig. |6] The CO 
with 1^2 (lower row) have bilayers with a net charge, indeed 
for (■iiO)/(000) the whole structure would be charged. 
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for each bilayer, the upper layer has a surplus of Fe'^^ ions (ji^y) [resp. (^i^O)] with the one from(00|) [resp. (000)] 
and the lower layer a surplus of Fe^+ ions, which implies 
an electric dipole moment between the layers (indicated by 
arrows), i.e. this CO is indeed ferroelectric, rendering the 
bilayers polar. 

4.1 Long-range CO from diffraction In the case of 
LuFe204, actual superstructure reflections were observed, 
below Tco ^ 320 K, by several groups using various 
diffraction techniques (Fig. |7]i. Most studies resolving 
3D peaks found the principal superstructure reflections 
to be near positions ( ^ , ^ ,ha lfinteger) and symmetry- 
equivalent |fT7l[T8l[T5]|fep3]|44l52 5354 555657 58 
l59l . with in-plane/out-of-plane correlation lengths given as 
e.g. 30/6 nm 153]] or 80/7 nm l(55l . Resonant x-ray diffrac- 
tion at the Fe K-edge of these reflections support their 
association with Fe'^+/^+ CO n7ll42ll54ll58l . Sharp re- 
flections at other positions, including near (i, i, integer) 
160116 Ill62l are sometimes observed in electron diffraction 
studies, generally on samples with higher oxygen content 
II59II62I . The additional reflections in 1591 were attributed 
to ordering of excess oxygen ions. 



(^^0)/(000) with Ti gives the ferroelectric CO with po- 
lar bilayers proposed in ||T7] . However, since the observed 
(also in ITTl ) superstructure reflections follow ( H |) prop- 
agation, one of the two CO shown on the right of Fig. 
|6] should be realized, either (top) the antiferroelectric CO 
with polar bilayers but zero net polarization, as proposed 
in l44l . or (bottom) the CO with stacking of oppositely 
charged bilayers. The latter was initially dismissed because 
it involves a transfer of charge between neighboring bilay- 
ers separated by ~ 6 A, but an experimental determination 
of which CO is realized requires a full structural refine- 
ment. 

This was recently achieved |T8l , the key being the 
screening of a large quantity of small crystals of the highest 
quality, as judged by M{T) (Fig. [3]l, for one with almost 
only one of the three possible CO-domains populated. The 
structure obtained at 210 K, readily refined in the space 
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Figure 7 Charge order CO superstructure in LuFe204 by different techniques, a) X-ray diffraction scattered intensity at 
(hhl) at 200 and 360 K. From g?), © 2008 American Physical Society, b) Neutron diffraction scans along (5^^) at three 
temperatures. From ll52l . © 2008 American Physical Society, c) Neutron diffraction intensity in spin-flip channel at (hhi) 
at 245 K. From (291, © 2012 American Physical Society, d) Electron-diffraction dark-field image taken at 92 K with a 
(■i i |) superstructure spot, circled in the inset. The CO domains extend ~ 30 nm in the a5-plane and ~ 6 nm along c. Inset: 
electron diffraction pattern with (110) incidence. From ll53l . © 2009 American Physical Society. 



group C2/m with a residual of 5.96%, contains four iron- 
sites with valences, determined with the bond-valence-sum 
(BVS) method 123, 1.9, 2.1, 2.8, and 2.9 HS). Given the 
uncertainties inherent to the empirical BVS-method, these 
values can be considered to be remarkably close to the 
ideal values of 2 and 3 (for a comparison of other charge 
ordering ferrites see Fig. 5 in |l68l). The thus identified CO 
is the one with charged bilayers, shown in Fig. |6]bottom 
right. Attempts to force-refine the other CO or to refine in a 
lower-symmetry space group confirmed this structure solu- 
tion, and independent confirmation was also obtained indi- 
rectly via spin-charge coupling (c.f. Sec. 15.3) . The ground- 
state CO of LuFe204 is thus not ferroelectric and does not 
contain the initially proposed polar bilayers. 

4.2 Short-range correlations and charge dynam- 
ics Above Tco, the sharp CO superstructure reflections in 
LuFe204 are replaced by diffuse scattering along {^-k^), 
with a characteristic zig-zag pattern, shown in Fig. |7^) 
lower graph, discernable up to about 550 K 14311611 . A 
detailed analysis of the scattered intensity (441 suggests 
that it can be described as very broad and overlapping 
peaks at positions of structural reflections ±(i — S',^ — 
S' , 0) and symmetry-equivalent. Corresponding to the large 




Figure 8 Example of a possible anti-phase-boundary 
(dashed) perpendicular to the in-plane propagation of CO, 
which switches the majority valence. Top-view of bilayer 
(•=Fe^+,o=Fe^+), Fe ions in the lower layer are drawn 
smaller. 



peak-width along i, at 360 K a correlation extending only 
to 2 — 3 bilayers was concluded, whereas the remaining rel- 
ative sharpness along (hhO) suggests that in each individ- 
ual bilayer, medium-range CO is maintained much above 
Tco- Applying symmetry-analysis on the positions of the 
diffuse peaks suggests a marginal tendency towards ferro- 
stacking of neighboring bilayers, one of the CO shown in 
Fig.|6]left. In (44l. representation Fi and thus a tendency 
towards ferroelectric CO, was assumed. However, as it is 
now clarified that in the long-range ordered state, repre- 
sentation ^2 with charged rather than polar bilayers is re- 
alized, it seems more likely that the short-range CO is also 
following Fi, without polar bilayers ever being present. 

This seems at first impossible, as this corresponds to 
an overall net charge of the structure. However, given the 
small correlation-volume deduced from the broadness of 
the diffuse peaks, a charge-transfer over a few bilayers 
(via a presently unknown pathway) is sufficient, and fea- 
sible considering that this is not a static, but a dynamic 
short-range CO, i.e. short-range-correlated valence fluc- 
tuations. Mossbauer spectroscopy can directly assess the 
dynamics of CO and charge fluctuations: spectra of sev- 
eral i?Fe204 were analyzed with a model featuring iso- 
mer shifts, quadrupole splittings, and linewidths for both 
valence states and additionally a Fe^^ -^ Fc"^^ hopping 
frequency, showing still discernable natural hopping even 
somewhat below Tco 125)1 . 



Fig. |9^) shows the hopping frequency iy\ 



loppm^ 



deter- 
mined in a study on powdered crystals of optimal quality. 
The hopping above Tco can be described by an Arrhe- 
nius law with an activation energy of ~ 0.16 cV. J^hopping 
decreases continuously upon cooling through Tco. and 
below can again be described by an Arrhenius law, but 
with a higher activation energy ~ 0.36 eV. The Arrhenius 
law with a change of slope at Tco is consistent with the 
T— dependence of the dc conductivity extracted by dielec- 
tric spectroscopy on similar samples ||T91 . The substantial 
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Figure 9 Charge dynamics in LuFe204. From 
1631 , © 2008 American Physical Society, a) 
Arrhenius-plot of the Fe^+ — >■ Fe''+ hopping 
frequency determined from Mossbauer spec- 
tra (o, bars are from KSI ). b) Temperature- 
dependence of the strength Weff of an optical ex- 
citation assigned to Fe^"*" — > Fe^^ charge trans- 
fer. 
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Figure 10 Electric-field effects in LuFe204. a) Current-voltage hysteresis loop at 300 K, with the current limited to 
100 mA. From [i69l . © 2008 American Institute of Physics, b) Resistivity p vs T-i/^, fuH symbols: low S-field d.c. 
measurements; open symbols: high-field pulsed measurements. From ITOl . © 2011 American Institute of Physics, c) 
Neutron diffraction intensity of ( |||) as a function of calibrated temperature, obtained by cooling with 0.7 mA (A) or no 
current (•). From 1561 . © 2010 American Physical Society. 



electron hopping present even below Tqq, where diffrac- 
tion indicates long-range CO, likely involves fluctuations 
of anti-phase boundaries, a scenario shown to be at work 
in another CO ferrite, and connected to T— dependent ap- 
parent incommensuration f65l . In LuFe204 the incommen- 
suration becomes constant only below ~' 170 K, suggesting 
that the charge dynamics become completely frozen only 
at temperatures substantially below Tco- 

The capacity of Fc^^ -^ Y(?^ charge transfer can be 
probed not only through natural electron hopping, but 
also by driving the process at optical frequencies. Fig. 
I9J5) shows the T— dependence of the strength of an op- 
tical excitation assigned to Fc'^^ -^ Fe'^^ charge transfer 



J. Upon cooling, it starts to decrease around 500 K, 
reasonably close to the reported I143II611 onset of diffuse 
scattering, and it bottoms out only at Tlt ~^ 170 K, consis- 
tent with the freezing of charge dynamics at this temper- 
ature. The nature of the phase below Tlt will be further 
discussed in Sec. 16.21 

4.3 Electric-field and current effects on CO 
Given that the pyroelectric current measurements shown 
in Fig. |5^) had been prepared by first cooling the sample 
in an electric field, and that diffraction analyses failed to 
deduce a ferroelectric CO, it had been proposed |l44l that 
such a ferroelectric CO might be stabilized by electric 
fields E or by cooling in E. However, at low T electric 
fields up to 20 kV/cm were shown to have no effect on 
the CO 1561 . At higher T, intriguing electric-field effects 
have been observed macroscopically in the form of non- 



linear current- voltage characteristics (Fig. [TOt ) 1169117 1 1 : 
Rather than eventually reaching a state with lower resis- 
tance as might be expected for a ferroelectric phase, a 
current-breakthrough is observed. The breakthrough de- 
pends on the environment so that a possible application 
in gas-sensing has been proposed |72j. Effects were also 
reported on dielectric constants 11731 and magnetism ||57] 
I74ll75l . With in-situ electron diffraction experiments, a 
disappearance of the CO superstructure reflections at this 
threshold was observed 07611771 . suggesting an electric- 
field or current-induced melting of the CO. In contrast, an 
ex-situ in-field cooling experiment suggests improved CO 
correlations, without affecting the intensity-distribution 
between different reflections Ii78j . 

However, current-voltage measurements using single 
short (^ 1 ms) current pulses llTOl found Ohmic behavior 
in polycrystalline LuFe204 (see Fig. [TOb ) up to E much 
higher than fields where the breakdown in dc measure- 
ments occurs, with the time-dependence on longer pulses 
suggesting self-heating effects fTOl . Similar pulsed mea- 
surements on YbFe204 fW\ suggest non-linearity in the 
in-plane conductivity, attributed to sliding charge density 
and possibly a break-down of CO, whereas large self- 
heating effects for longer pulses are also found. An in-situ 
neutron diffraction experiment ll56l demonstrates the ab- 
sence of any current-dependence on either charge or mag- 
netic order when the sample-temperature is properly re- 
calibrated (Fig. [Tot). A CO unaffected by static (melting 
of the CO by 1.55 eV photons was demonstrated by a re- 
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cent pump-probe experiment lISOl ) electric fields or cur- 
rents was also found by in-situ infrared spectroscopy llSTI . 
by in-situ x-ray diffraction with the real sample tempera- 
ture established from the thermal expansion of the silver 
electrodes 1821 . and by resonant x-ray diffraction after in- 
field cooling Il58l . Although a recent infrared-spectroscopy 
study |83| suggests very subtle field-effects not attributable 
to self-heating, the conclusion that a different, possibly fer- 
roelectric, CO cannot be stabilized by electric fields is un- 
avoidable. 

4.4 Theoretical considerations The experimen- 
tally deduced CO with charged bilayers is surprising, be- 
cause CO is usually assumed to be driven by the repulsion 
of electrons between different sites ("Wigner crystalliza- 
tion"), which should result in minimizing the occurrence 
of neighbors of the same valence. It is clear that packing 
electrons close together in a charged bilayer should lead to 
a significant energy penalty. Although the full description 
of the 3D CO requires at least the four interactions Ui — C/4 
indicated in Fig. [T];) 01511161 . the clear indications by the 
diffuse scattering above Tco (Sec. 14.21 ) that correlations 
within the bilayers are much stronger than those between 
different bilayers make a first analysis based on only Ui 
and U2 (intrabilayer interactions) useful. Concerning the 
intralayer interaction Ui all the CO of Fig. |6] as well as 
many others, such as a stripe-like CO, are exactly degen- 
erate due to geometrical frustration (c.f. Sec. |2). However, 
concerning the three nearest neighbors (NN), which are 
in the other layer {U2), the charged bilayer CO (Fig. |6] 
bottom) have on average | same-valence NN, whereas 
the polar CO (Fig. |6] top) have only | same-valence NN, 
and is therefore more stable for repulsive 1/2- However, a 
stripe-like CO can have 1 same-valence NN, which should 
make it even more favorable, as found also by Monte-Carlo 
simulations taking into account the long-range nature of 
the Coulomb interaction ll84l . This is an indication that 
Coulomb repulsion is not the sole relevant driving force 
for CO in i?Fe204. It is not really surprising; a proof of 
purely electrostatically driven CO has been elusive in bulk 
materials and strong influence of lattice distortions on CO 
are well known, e.g. in manganites ll85l . In addition to 
lattice effects, indicated to be relevant in LuFe204 by in- 
frared ||7]|82 and Raman |83 spectroscopy, CO driven by 
magnetic exchange has also been demonstrated | 



Density-functional-theory (DFT) calculations per- 
formed for LuFe204 1841 suggest a lattice-contribution 
of about 40% to the total energy gain by CO for the fer- 
roelectric (402mcV/f.u.) and stripe (384mcV/f.u.) CO 
investigated, estimated by separately stabilizing the CO 
with fixed atomic positions and with relaxing the structure. 
Charged bilayer or antiferroelectric CO could not be ob- 
tained, because a cell with 3 bilayers was used. A ferrielec- 
tric CO with a fit stacking of bilayer polarizations (and a 
net polarization close to the one indicated in Fig.|5^, lend- 
ing support to this model) and later using a larger cell an 
antiferroelectric CO ll44l was found to be yet more stable, 



by a total of 15 mcV/f .u.. These small energy differences 
between antiferroelectric, ferroelectric, and stripe CO are 
an indication of the relevance of geometrical frustration, 
as expected (Sec. |2]). Spin-orbit coupling and magnetism 
was considered in a later study ll89l (see Sec. 15.3b . In this 
study, a fixed ferroelectric CO was used, but the reported 
energy gain of 78 mcV/f .u. is considerably larger than the 
above energy-differences between different CO, indicating 
that magnetism may be relevant in stabilizing the CO. The 
charged bilayer CO has yet to be tested by DFT. 

The CO in i?Fe204 has also been studied by lattice 
gas models. Modeling a single bilayer I90ll9ni92ll93ll94l 
already yields several competing CO phases, among them 
the ferroelectric CO (Fig. |6] top left) and non-polar CO 
phases with in-plane propagations such as (■^0)A\j)A^^) 
and (^^), again pointing out the importance of geomet- 
rical frustration. Because in |90.;9lJ92,!93;i94j charge neu- 
trality of this bilayer was implied, the charged bilayer CO 
was not considered. A lattice gas model defined on the 
whole 3D lattice of Fe has been considered in I15II16I . 
though without considering in detail the charged bilayer 
CO found experimentally. Starting point is the Hamilto- 
nian 



where i/j enumerate the primitive unit cells, i^/u' = 1,2 
denote the Fe ion in the primitive cell, a assumes the value 
1 (—1) for the valence state 2+ (3+) of the Fe ion, and 
U denotes screened Coulomb-interactions between differ- 
ent Fe ions. In reciprocal space, the eigenvalues of the 
Fourier-transformed Hamiltonian dl) give the energies for 
charge patterns defined by a corresponding eigenvector and 
a propagation vector, and the charge pattern with the lowest 
energy corresponds to the realized CO. Four interactions 
(Ui — f/4 in Fig.[T];) are considered. It is clear that at least 
interactions Ui — C/3 are needed to yield a 3D CO: They 
lead to instabilities at (hhO^ but these are degenerate along 
£, yielding only 2D CO 1151 . This is due to the rhombohe- 
dral stacking of the bilayers |[T6l . a geometrical frustration 
effect. 

With a non-zero C/4, depending on the signs of C/4 and 
C/2 ■ C/3 possible ordering wavevectors (^ ±S,^ ±6, |) or 
(i ± (5, i ± (5, 0) are obtained. C/4 < (attractive, i.e. favor- 
ing same valence) must be assumed to yield the observed 
propagation vector, which was proposed to be due to over- 
screening IllSj . The observed (c.f. Sec. 14.2b discrepancy 
between the dominant fluctuations above Tco with peaks 
near (i, i, integer) and the long-range CO below Tco with 
reflections near (i, i, half-integer) is difficult to reconcile 
with a mean-field lattice gas model. In lT6l . the coupling to 
a non-critical mode, specifically a particular Raman-active 
phonon-mode, was proposed. Raman scattering has so far 
been conducted only on polycrystalline samples, not find- 
ing any indications for this scenario Il87l . However, a prob- 
lem of the model is that with the simple quadratic Hamilto- 
nian ([T]), an instability at only a single q-position can be ob- 
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tained, i.e. the additionally observed reflections near (00|) 
cannot be explained. 

4.5 Ion-size effects on CO Given the large distor- 
tions of the bilayer in i?Fe204 as a function of R ion- 
size (Sec. 12] and Fig. [T}?), different CO for different R 
may be expected. Yb is quite close in ion size to Lu, and 
whereas long-range CO was not yet found, presumably due 
to stoichiometry -problems, diffuse scattering near (i-i£) is 
generally observed [95 1, and an electron-diffraction study 
1961 even found maxima at £ half-integer. Given similar 
spin-charge coupling (see Sec. l5.3l l, the same CO as for Lu 
is likely. For R=ln, which has a significantly smaller ion 
size, single crystals are not available, but synchrotron x-ray 
powder diffraction suggests the appearance of long-range 
CO below ~ 250 K, with peaks indexed by (H^), with 
r incommensurate fTT): the limit of thick bilayers seem 
to favor an incommensuration out-of-plane rather than in- 
plane. 

For the other extreme of thin bilayers, R =Y, the same 
type of diffuse scattering along (-^^t) is observed at all 
T (e.g. 1971 ) for samples with significant oxygen defi- 
ciency as in off-stoichiometric Lu or Yb samples. How- 
ever, the behavior of sufficiently stoichiometric material is 
more complex. Abrupt changes of resistivity accompany- 
ing the two magnetic transitions marked by arrows in Fig.|2] 
l98l suggest these are also structural transitions, with mon- 
oclinic and finally triclinic distortions as indicated by pow- 
der x-ray diffraction [,99J . 

Because weak superstructure peaks are difficult to ob- 
serve in powder diffraction, and sufficiently stoichiomet- 
ric single crystals were absent until recently, subsequent 
work on the CO of stoichiometric YFe204 focused on 
electron-diffraction (e.g. jgTlfTOOlfToTl ). At low T, these 
studies found clear superstructure spots with very differ- 
ent propagation vectors, such as (-^|-L.) between 85 and 
130 K, and consistent with triclinic symmetry. At higher T 
often a (iiO)-type superstructure is observed, and some- 
times (5^0)- or more "conventional" (ii|)-types, partly 
with different phases coexisting or relaxing from one to the 
other as a function of time. The question of 3D CO in sto- 
ichiometric YFe204 above the magnetic ordering temper- 
ature, in particular at room temperature, is not clear, with 
either only diffuse {■^^i) lines reported II 1001 or with coex- 
istence of diffuse lines with broad (hk^) spots, the latter 
by dark-field imaging assessed to correspond to 3D CO do- 
mains correlated over some 30 bilayers ll97l . 

These observations suggest a subtle competition of sev- 
eral CO phases, that correspond to highly complex charge 
patterns, in YFe204, reminiscent of some of the predic- 
tions based on single-bilayer model Hamiltonians (see Sec. 
I4.4l i. To elucidate the strong difference in CO to LuFe204 
will require detailed modeling, but first the real-space 
CO giving rise to the observed superstructure reflections 
will have to be solved, presumably by single-crystal x-ray 
diffraction on sufficiently stoichiometric crystals II3TI. 
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Figure 11 Magnetic 
field H dependence 
of magnetization M 
and the intensity of the 
(T02)+(00|) magnetic 



reflection 
LuFe204. 



at 220 K in 
After l29l. 



4.6 Fe-site substitution and intercalation As 

mentioned in Sec. |2l possible changes in oxygen sto- 
ichiometry hamper substitution studies on the Fe site. 
Nevertheless, the expected quick depression of CO or 
charge diffuse scattering seems to hold (see, e.g., II102I ). 
For intercalated materials, e.g. Lu2Fe307, superstructure 
reflections have been observed by electron diffraction, in- 
terestingly up to 750 K, and interpreted in terms of both 
the bilayers and the intercalated single Fe-layers exhibiting 
CO O8lll03l[T04l . X-ray and neutron diffraction on single 
crystals would certainly be desirable. 

5 Magnetism and Spin-Charge coupling The 

i?Fe204 family first attracted attention due to intriguing 
magnetic properties of Ising spins (pointing perpendicular 
to the layers) on a strongly frustrated lattice, and the mag- 
netic properties depend strongly on oxygen-stoichiometry 
(c.f. Sec.HI. 

Off-stoichiometric samples typically show signifi- 
cantly higher magnetization when cooled in a magnetic 
field than when cooled in zero field (c.f. Figs. |2]and[3t), 
with a remanent magnetic moment after cooling in strong 
magnetic fields between 1.5 (Y) and 2.8/iB/f.u. (Lu) 
11051 . Frequency-dependence in ac magnetic susceptibil- 
ity is observed below about 220 — 250 K, with detailed 
analyses suggesting either spin-glass- 01061 or cluster- 
glass-like [107,108,1091 freezing, with large clusters of 
size ^ 100 nm in-plane directly observed by magnetic- 
force-microscopy in one example 01071 . Neutron diffrac- 
tion studies suggest 2D magnetic correlations, with diffuse 
magnetic scattering typically observed along the {-kkO 
line 130111 101 similar to the CO in non-stoichiometric sam- 
ples (c.f. Sec.|4|i. This suggests that mainly the correlation 
between different bilayers is lost. In less off-stoichiometric 
samples, a tendency towards a development of peaks along 
this line has been observed l30l . but so far few neutron 
diffraction studies of stoichiometric crystals has been 
performed, and to date only the magnetic structure of 
LuFe204 is known, which will therefore be the focus in 
the following. 

5.1 3D magnetic correlations in LuFe204 Sharp 
magnetic Bragg-reflections have been observed after cool- 
ing below Tn ^ 240 K at (i, i,integer) and (i, i,half- 
integer) positions in several single crystals of LuFe204 by 
neutron diffraction 029)15 2115 5 115 6ll57l and soft resonant x- 
ray diffraction ll29l. Magnetic correlations between differ- 
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Figure ^2 H - T Magnetic 
phase diagram of LuFe204 from 
data of I29II63I . Shown are a 
low-_ff antifeiTomagnetic (AFM) 
and a high-iif ferrimagnetic (fM) 
phase near the ordering temper- 
ature Tn ~ 240 K, as well 
as a lower-temperature phase 
(LT) with re-entrant disorder 0521 
and a subtle structural modifi- 
cation 1631 . Hysteretic regions 
where more than one phase can 
be stabilized are hatched. The 
LT -)■ fM transition line (A) 
reaches 7 T at 50 K (not shown) 
1631 . The spin structures of the 
AFM and fM phases (291 are 
also shown, with large semi- 
transparent arrows indicating the 
bilayer net magnetizations. 



ent bilayers have also been deduced from powder neutron 
diffraction 11107111111 . An example of (^^i) scans at dif- 
ferent temperatures is shown in Fig. |7j)): the weak reflec- 
tions at (i i , half-integer) at 280 K (thick olive) are due 
to CO, the much larger intensity mainly at ( -i , -i ,integer) at 
220 K is magnetic in origin, as shown by polarized neutron 
diffraction (see, e.g. Fig. 5 in ||29l ), with an out-of-plane 
correlation length of ISmii. At 130 K (red) a large dif- 
fuse background appears together with peak-broadening, 
suggesting a re-entrant magnetic disordering II52II561 . This 
low-temperature state (LT phase) will be discussed in Sec. 
6:21 



In contrast to stoichiometric YFe204, where the ab- 
sence of any magnetic-field induced phase transitions is 
reported II27II991 . magnetization-measurements on sto- 
ichiometric LuFe204 reveal a complex magnetic field 
(i?||(001)hex)-temperature (T) phase diagram 112911631 . 
shown in Fig. [12] Apart from the already mentioned LT 
phase with reentrant magnetic disorder, and the param- 
agnetic state, two further phases are present. Isothermal 
magnetization M{H) loops (e.g. Fig. [TT] grey line) are 
consistent with a hysteretic metamagnetic transition be- 
tween a low-H phase that is antiferromagnetic (AFM) 
and a high-iJ phase that is ferrimagnetic (fM), with a net 
magnetic moment at low temperatures (the fM phase re- 
mains stable down to T = when cooled in hqH > 2 T) 
of ^ S/is/f.u., suggesting a fti-like arrangement of Fe 
spins. Partially strong hysteresis of the transition between 
AFM and fM phases (hatched in Fig. [T2] |. also present 
between LT and either fM or AFM, is one of the unusual 
features of the phase diagram. Magnetic Bragg reflections 
are observed for both the AFM and the fM phase at i^-^i), 
with an intensity shift surprisingly from i integer to half- 
integer reflections going from AFM to fM 112911551 . As 
expected from the presence of a net moment, in the fM 
phase magnetic intensity is also observed on structural 



reflections. In contrast, magnetic intensity at (00|)-type 
reflections is a feature of the AFM phase (see Fig.flTTi. 

5.2 Competing spinstructures Given the corre- 
spondence of Ising spins t or | with valence states 2+ 
and 3+, both being binary orders, the most straightforward 
way to analyze spinstructures is to neglect the different 
spin quantum numbers for the different valence states of 
the Fe ions and proceed from the hexagonal cell analogous 
to the CO. As the same propagation vectors are involved 
as considered for the CO in Sec. 14.1! there will again be 
three domains and the same four possible orderings shown 
in Fig.|6]also then apply to the spin order (SO), with black 
and white coloring now indicating spins t and |. If the 
weaker (-i, -i , half-integer) magnetic reflections are consid- 
ered to be due to a decoration by the CO, then the proper 
SO should be the ordering shown in Fig. |6] bottom left, 
because for a Fi representation there would be no inten- 
sity on (HO), which is the strongest magnetic reflection 
(Fig. I7J5). This ferrimagnetic SO was found to fit well 
all (i, ■i,integer)-type reflections and was thus initially 
proposed to apply ll52l . 

However, the later found ||29| magnetic intensity at 
(10|) (Fig.fTTTi is inconsistent with this SO, showing that 
symmetry-analysis based on the hexagonal cell cannot de- 
scribe the SO of LuFe204. In retrospective this is not sur- 
prising, because the SO takes place at Tn < Tco and if any 
spin-charge coupling is present (c.f. Sec. 15.3) the proper 
cell on which SO is to be considered is the CO supercell. 
Because the magnetic cell is equal to the CO cell, sym- 
metry analysis is then based on a (000) propagation, and 
with four independent sets of Fe sites present in this cell, 
and four irreducible representations for each site, many 
SO have to be considered. According to Landau-theory, all 
sites should order according to the same irreducible rep- 
resentation if the ordering occurs in a second-order phase 
transition at a single temperature. This condition yields 40 
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possible SO, which however were all found to be inconsis- 
tent with the experimental data 1291 . In a brute-force ap- 
proach all 3^2 = 531'441 SO of Ising spins (t,|,0, al- 
lowing for some Fe sites remaining disordered) were com- 
pared with the data for both AFM and fM phases, which 
identified the SO shown in Fig.[T2]as the correct ones, ver- 
ified by refinement |29| . For these SO, all four sets of Fe 
sites order according to an irreducible representation, but 
with different ones for different sets, implying that the SO 
does not occur via a single second-order phase transition. 
Indeed, given that both AFM and fM, separated by a first- 
order transition, are present near Tn and H = Q, a simple 
second-order transition seems unlikely. 

AFM and fM SO are closely related. In both cases, each 
bilayer (and each individual layer) has a 2 : 1 (tt^ or H.'f) 
ferrimagnetic arrangement of the Fe spins, with the bilayer 
net moment indicated by thick semi-transparent arrows in 
Fig. [121 For the ferrimagnetic fM phase, the net moments 
of the two bilayers of the cell are both '[, yielding a net 
moment, as observed. For the AFM phase, the spins in one 
of the two bilayers are exactly as in the fM phase, and for 
the other bilayer they are all reversed, yielding the antifer- 
romagnetic behavior indicated by magnetization measure- 
ments. Intriguingly, as can be seen in Fig. [12] the AFM-fM 
transition extrapolates to iJ « for T -^ Tn, implying 
the near-degeneracy of these two SO phases at Tn in zero 
field, a hallmark of geometrical frustration. A competition 
of AFM and fM SO has also been observed in a study on 
polycrystalline samples [1111|. Because AFM and fM SO 
are distinct only by the flip of all spins in one of the bilay- 
ers, this suggests that the spin correlations within the bilay- 
ers are much stronger than the inter-bilayer correlations. In 
analogy to the CO (Sec. l4.2l i. a random stacking of the net 
moments of still relatively well spin-ordered bilayers can 
then be expected above Tn, with diffuse magnetic scat- 
tering along (55^). This has indeed been observed (Fig. 
[7};), up to room temperature, and deviations of the mag- 
netic susceptibility from Curie-Weiss behavior suggest that 
spin-fluctuations correlated within each bilayer remain rel- 
evant beyond even Tco lf29| . Note that these correlations 
are dynamic (faster than ^MHz) rather than static, because 
no indications of magnetic order are visible in Mossbauer 
spectra at 260 K 1631. 

5.3 Spin-charge coupling 

Experimental indications Experimental indica- 
tions of a spin-charge coupling by anomalies in the T- 
dependence of the intensity or incommensuration of the 
CO superstructure reflections at T^ have been observed 
I44II58I . but they are weak. However, this might be ex- 
pected if the SO within individual bilayers is formed al- 
ready above T^. The SO refinements f29\ yielded some 
hints of a strict assignment of SO and CO, but the different 
magnetic moments of Fe^+ and Fe'^"'" ions could not be 
reliably determined. However, spectroscopic techniques 
performed on a variety of samples provide strong evidence 
that such spin-charge coupling is a common feature at 
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Figure 13 X-ray absorption spectra XAS with right (/^+) 
and left (fi-) circular polarized photons and X-ray Mag- 
netic Circular Dichroism XMCD Aji = ji^ — ji^ of 
LuFe204 across the Fe £2/3 edges in 4 T (perpendicular 
to the layers) at 120 K and 260 K (inset). Also shown is the 
integration of the XMCD, with the sum rules yielding an 
orbital magnetic moment of ^ 0.7 /^s/f .u.. From lITSl . @ 
2012 American Physical Society. 



least in YbFe204 and LuFe204 l[T8ll64l [n2l[Tl3llll4l . In 
the first studies, polarized Mossbauer spectroscopy was 
applied to single crystals of YbFe204 |113j and LuFe204 
[114:1 under varying H, it was concluded that all of the 
Fe^"*" spins as well as i of the Fe^+ spins were aligned 
\\H, with the remaining | of the Fe^+ spins aligned an- 
tiparallel to H. 

Later studies by x-ray magnetic circular dichroism 
(XMCD) confirmed this arrangement for LuFe204 ITSl 
I64llll2l . Fig. [13] shows the absorption spectra and XMCD 
measured (in the fM-phase) on a crystal of the same type 
as used to establish the magnetic phase diagram and SO 
(c.f. Fig. \\% . The larger negative peak of the XMCD sig- 
nal at the Fe^+ L3 edge and the smaller positive peak at 
the Fe'^^ L3 edge directly implies a larger Fe^"*" net mo- 
ment in magnetic field-direction, and a smaller Fe'^^ net 
moment opposite to this, which given the known magnetic 
cell-size and saturation magnetization leaves the above 
spin configuration as the only possible one. From the 
XMCD sum rules, a sizeable orbital magnetic moment of 
- 0.7 - 0.8Aii3/Fe^+ ion can be concluded l29ll64l[TT2l . 
rendering the Fe^+ total magnetic moment to ~ 4.75 /is, 
quite close to the spin-only Fe'^+ moment of 5 jib- The re- 
sulting saturation moment of 4.75 — -^ • 5 ~ 3.08 /is/f .u. 
is in agreement with magnetization data. 

In Fig. [12] both spin and charge structures are shown, 
and it can be readily seen that for the fM-phase the align- 
ment of Fe-^+ and Fe'^^ moments according to XMCD and 
Mossbauer spectroscopy is indeed observed, i.e. the previ- 
ously discussed CO and SO are consistent with strict spin- 
charge coupling as deduced spectroscopically. Moreover, 
the SO shown in Fig. [T2l combined with the spectroscopic 
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Figure 14 Intensity of the (|||) CO superstructure peak 
measured by neutron diffraction at 300 K, after cooling in 
magnetic fields of different strengths. From [35], © 2009 
American Physical Society. 



results is incompatible with any of the CO shown in Fig. 
|6]except the charged bilayer stacking deduced in Sec. 14. II 
providing additional support for this model of CO IITSl . 

The inset of Fig. [T3] shows a XMCD spectrum mea- 
sured in 4 T above Tn, at 260 K. Though much reduced in 
amplitude, the shape is still the same as in the fM-phase 
below Tn. This is consistent with the model of randomly- 
stacked but otherwise still well spin-ordered bilayers in the 
paramagnetic state: Polarized by a sufficiently strong mag- 
netic field, the random-stacking of 50% t and 50% I bi- 
layer net moments will be broken, with the overall spin 
arrangement tending towards the fM SO. It also shows that 
the spin-charge coupling is very robust, not requiring long- 
range order The magnetic correlations could thus be im- 
portant not only for the magnetic transition at Tn, but also 
for the charge ordering at Tco- 

Interestingly, a magnetic-field influence on the CO has 
been observed for LuFe204 1551 : Fig.[T4lshows the inten- 
sity of the (| ||) CO superstructure reflection, measured at 
300 K, after cooling the crystal in different magnetic fields 
from 350 K. A significant suppression of ~ 15% already 
by rather small fields is visible. This is surprising, given 
that at the AFM-fM metamagnetic transition no change of 
the CO was observed L18J . However, in contrast to these 
measurements with H\\c, in lf55l the field was applied in 
[110] direction. This in-plane H breaks the three-fold rota- 
tion symmetry of the crystal, and likely affects the relative 
stability of the three magnetic domains (based, as the CO 
domains, on the three monoclinic cells shown in Fig. [T^). 
The strict spin-charge coupling implies the correspondence 
of magnetic and CO domains, and therefore an in-plane H 
should affect the domain population. The (fH) reflection 
studied in ll55l belongs to one of the CO domains, its in- 
tensity suppression by H suggests that population of other 
CO domains is favored, which could be tested by simulta- 
neously tracking different CO reflections belonging to all 
domains. 

Theoretical considerations The single-bilayer lat- 
tice gas models Il90ll9ni92ll93ll94l also consider differ- 



ent SO and spin-charge coupling, finding that the spin- 
correlations tend to stabilize CO with (i-i^)-type propa- 



gation compared to alternatives (c.f. Sec. I4.4l i. However, 
because these models do not consider the charged bilay- 
ers indicated by experiment, the focus below is on more 
phenomenological considerations and the DFT calculation 



Because the inter-bilayer spin-coupling is very weak 
(c.f. Sec. I5.2l l. it is a good first approximation to con- 
sider each bilayer, the one with Fe^+ majority and the 
one with Fe'^^ majority, individually. Various magnetic su- 
perexchange parameters between Fe ions of the same layer 
and between the two layers have then to be examined, 
which may be expected to strongly depend on the valence 
state of the respective Fe ions. An experimental determi- 
nation of exchange parameters could be done by inelastic 
neutron scattering, which so far have been reported only in 
a limited range on a crystal [78 1, and on a polycrystalline 
sample, with the observation of a spin-gap ^ 7 mcV con- 
firming the Ising nature of spins 01111 . 

The exchange parameters were estimated theoretically 
by mapping the energy differences between differently or- 
dered states obtained by DFT to the corresponding en- 
ergy differences of an Ising Hamiltonian ||89l . This study 
considered an individual bilayer, however, one assumed 
to exhibit the ferroelectric CO (Fig. |6] top left). The by 
far largest interaction was found to be an antiferromag- 
netic superexchange between Fe'^+ ions of the same Fe'^^ 
majority layer From this, an in-plane antiferromagnetic 
arrangement of Fe'^^ moments is expected, which is in- 
deed observed experimentally (see Fig. [12). A correspond- 
ing calculation taking into account the charged bilayer CO 
would be highly desirable: as mentioned already in Sec. 
14.41 the large energy-differences between different spin- 
states found in |f89]| imply that the CO and SO have to be 
calculated together, in a cell of sufficient size, which makes 
the calculations challenging. 

5.4 Off-stoichiometric samples and other rare 
earths The agreement between XMCD experiments con- 
ducted on samples with f l8j and without ||64| long-range 
CO indicates that the strict spin-charge coupling with the 

Pg2+ ^^^ /'pe^+ -f^^ arrangement is a robust and generic 
feature of LuFe204. Because the superexchange interac- 
tions leading to spin-charge coupling and SO depend on 
the underlying CO, it seems very likely that in all samples 
individual bilayers exhibit both CO and SO as discussed 
in Sees. 14.11 14.21 and 15.21 with only the occurrence of 
long-range 3D order with correlations between different 
bilayers being affected by stoichiometry. The generic ob- 
servation of at least diffuse charge (e.g. 01071 ) and diffuse 
magnetic (e.g. fTTOl) scattering at the (kkO line supports 
this assessment. Because the charge dynamics slows down 
upon cooling regardless of stoichiometry ll25l . a glassy 
freezing of the essentially still charge-ordered bilayers 
can be expected to replace the long-range CO when the 
latter is prevented. Similarly, a glassy freezing of SO bi- 



Copyright line will be provided by the publisher 



pss header will be provided by the publisher 



13 



layers can be expected to replace the long-range AFM 
phase upon cooling, and such a cluster-glass-like behav- 
ior involving frequency-dependence in ac susceptibility 
has indeed been found for some samples, e.g. in 01081 . 
On the other hand, when a magnetic field forces the net 
moments of all bilayers to align, a stable fM-like phase 
with the same saturation moment should be recovered, as 
is indeed suggested by all macroscopic measurements (see 
e-g- ©61 1081 for corresponding phase diagrams) and by 
electron holography 01151 . Note, however, that because of 
the underlying random stacking of Fe^+ and Fe'^+ majority 
bilayers, this fM phase is not expected to be 3D ordered, 
but should rather exhibit diffuse magnetic scattering along 
(■||^), which could be checked by in-field neutron diffrac- 
tion on less-stoichiometric crystals. Thus, based on the 
above, for off-stoichiometric samples, exactly the "spin- 
glass-like 2D-ferrimagnetic order" proposed early on [11101 
is expected to occur 

The above considerations should hold also for YbFe204, 
which has a similar R ion size. Indeed, not only the same 
diffuse scattering at (■^■^i) has been observed 0951 as on 
non-stoichiometric LuFe204, but the same spin-charge 
coupling is also deduced I113L For i?=In with consider- 
ably smaller ion size, magnetism and spin-charge coupling 
appear to be rather different as indicated both by a much 
smaller saturation moment of l.S/^s/f.u. and dissimilar 
Mossbauer spectra [[ijj. An XMCD study on InFe204 
might be helpful to elucidate spin-charge coupling in this 
material. For R with considerably larger ionic radii than 
Lu or Yb, there are also clear differences. For example, 
in stoichiometric YFe204 there is no fM phase, and given 
the wildly varying CO observed at low T (see Sec. I4.6l l 
a similar spin-charge coupling as in LuFe204 is rather 
doubtful, despite similarities of diffuse scattering observed 
in off-stoichiometric crystals. In agreement with substan- 
tial changes upon increasing the rare earth ionic size, an 
XMCD study on LuFe204 with Lu partially substituted 
by the larger Y or Er ions 111161 , found marked changes in 
the spectra. Whereas the negative peak at the Fe^+ edge 
remained largely unaffected by the substitution, suggest- 
ing that "all Fe'^^ spins aligned with i/" still holds, the 
positive peak at the Fe'^^ edge was completely suppressed 
for ^ 50% substitution, suggesting an antiferromagnetic 
order of Fe^^^ spins. 

6 The ground state: orbital order? 

6.1 Orbital order, orbital magnetic momentum, 
and magnetic anisotropy Concomitant with CO and 
SO, orbital order (OO) typically occurs in other transition 
metal oxides, strongly influencing the magnetism as well 
II117I . The possible presence of OO has been considered 
for LuFe204 theoretically [|84„;90,|92| , for example the 
DFT calculations 1841 suggest a clear OO. However, the 
discrepancy of the predicted CO with the experimentally 
deduced one suggests that the predictions of OO also have 
to be taken with care. 
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Figure 15 Sketch of the energy-splitting of the Fe 3d or- 
bitals (right) due to the trigonal-bipyramidal coordination 
(left). For Fe'^^, each orbital is singly occupied, the extra- 
electron of Fe^+ can go to either of two e" orbitals. 



The natural starting point for the consideration of OO is 
the trigonal bipyramidal crystal-field level diagram shown 
in Fig. [15] ForFe'^^, all five orbitals are singly occupied by 
a spin t electron, and there is no orbital degree of freedom. 
However, the extra ^-electron at Fe^"*" ions can go to an 
orbital of the doubly-degenerate lowest level. It should be 
mentioned that e' and e" levels are quite close in energy 
in LuFe204: while band-structure calculations suggested 
that e' is the lowest level 11841 . polarized x-ray absorption 
spectroscopy at the O if— edge has later shown that the e" 
level is actually lower Il64l . 

OO would then imply a preferential occupation of ei- 
ther the yz or the zx orbital — or a given real linear com- 
bination of these — for each site. An alternative to this 
is the preferential occupation of a complex linear com- 
bination that could carry an orbital angular momentum 
L ((-01 L |-0) = for tp real) and thus an orbital mag- 
netic moment. From a linear combination of yz and zx 
orbitals, no eigenstates of L^ or Ly can be obtained, but 

\/ 1 ("F-zs; — i ■ yz) is an eigenstate of Lz with eigenvalue 
771 = ±1, allowing an orbital magnetic moment of up to 
±l/iB/Fc^^ ion. The fact that an orbital magnetic moment 
~ 0.7 — 0.8 /.ts/Fc^^ is observed (see Sec. 15. 3b implies that 
in the fM state (and at low T for samples without long- 
range SO) OO cannot occur Consistent with this, no addi- 
tional reflections have ever been found, and no anisotropics 
in resonant diffraction suggestive of ferro-type OO have 
been found ll42l . For samples with long-range SO, this con- 
clusion can also be extended to the AFM and paramag- 
netic phases (c.f. Fig. fT2l i. because high-resolution x-ray 
diffraction indicates no structural changes that necessar- 
ily would accompany an OO-transition at the fM-AFM or 
AFM-paramagnetic transitions OTSl (though for the latter 
small intensity anomalies are suggested in II44II58I ). 

From the presence of an orbital magnetic moment nec- 
essarily perpendicular to the layers, the large magnetic 
anisotropy making this a model Ising-system follows im- 
mediately, because spin-orbit coupling aligns the direc- 
tions of spin- and orbital magnetic moments 11641 . The large 
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magnetic anisotropy is to a good part responsible for the 
giant coercivity of LuFe204, with coercive fields reaching 
around 10 T at low T 123111071 . and to "exchange-bias"- 
like properties 1 1 1 81 . 

6.2 The low-temperature phase The only phase in 
the H — T phase diagram (Fig.[T2]i consistent with OO is 
the LT-phase, because both AFM o LT and fM O LT 
transitions are also structural transitions as indicated by 
Bragg-peak splitting ll63l . As discussed in Sec. 14.21 the 
transition to this phase is accompanied by the complete 
freezing of residual charge dynamics, which is consistent 
with an OO scenario, given that OO often accompanies 
metal-insulator transitions 111 171 . Due to the impact of OO 
on the magnetic exchange interactions, a change of the 
magnetic structure upon entering the LT phase would then 
be expected, consistent with the observed re-entrant disor- 
der (Fig. Eh). 

These indications that OO might be the driving force 
behind the transition into the LT phase are, however, at 
most tentative. As discussed in Sec. 16.11 a full OO would 
imply the quenching of the orbital magnetic moment and 
the removal of the strong tendency of spins to be aligned 
perpendicular to the layers. The relatively good description 
of the ^—dependence of the diffuse component of magnetic 
scattering at 130 K assuming spin direction perpendicular 
to the layers Il52l is therefore an argument against an OO 
scenario. 

Structurally, despite of the Bragg-peak splitting, there 
is no basic change of the CO across Tlt: the interbilayer 
correlations of CO seem to improve slightly 1.44 J . but the 
cell size (no extra reflections), C2/m symmetry, and the 
valence states of the Fe ions estimated by bond-valence- 
sum are unaffected within error bars IITSl . Although small 
changes in lattice parameters have been reported, the 
changes of refined atom positions are much too subtle 
to conclude about a possible OO IfTSl . Further structural 
refinement studies across Tlt are certainly desirable. A 
further possibility to detect OO is to assess asymmetries in 
resonant x-ray diffraction, ideally at the Fe L edges at the 
(rT|) reflections 082111191 . although the analysis is com- 
plicated by the possible emergence of magnetic intensity 
at these positions in the LT phase 152115611 . 

The nature of the LT phase, and whether it is connected 
with OO therefore is an open question left for future work. 
Its experimental investigation is also hampered by its very 
fragile nature, as it is present in some samples showing 
long-range CO and SO II 1 8ll29ll52ll56ll63l (indications are 
also visible in M{T) on the polycrystalline samples stud- 
ied in 122111111 ), but not in others II15II55I . Due to this, it 
is also not completely clear whether the LT phase repre- 
sents the true ground state of LuFe204, though it appears 
likely. If the LT phase is the true ground state, it will make 
theoretical investigations with zero-temperature methods, 
such as DFT, difficult because of the presence of signifi- 
cant magnetic disorder 



7 Summary and outlook In the foUowing, I wifl 
summarize the main results reviewed and highlight critical 
points for future investigations on rare earth ferrites, con- 
cluding by answering the title question and by looking at 
the larger picture of possible CO-based ferroelectricity in 
other compounds. 

For the most studied rare earth ferrite, LuFe204, a con- 
sistent and rather comprehensive picture of strongly cou- 
pled charge- and spinorder has emerged from experimen- 
tal work. Single-crystal x-ray and neutron diffraction re- 
vealed sharp charge order (CO) superstructure peaks below 
Tco ~ 320 K and spin order (SO) superstructure peaks be- 
low Tn ^ 240 K in stoichiometric samples. Both are unaf- 
fected by static electric fields or currents, except for self- 
heating effects. The corresponding orders have been de- 
termined by refinements. The CO involves charged, rather 
than polar, bilayers. As revealed by Mossbauer and optical 
spectroscopies, substantial charge fluctuations occur even 
below Tqo, down to Tlt ~ 170 K. Between Tlt and Tn 
two related SO phases compete, due to geometrical frus- 
tration. Below TifT, a subtle structural distortion and re- 
entrant magnetic disorder occur Whether this is connected 
with possible orbital ordering is a major question left for 
future work. Mossbauer spectroscopy and x-ray magnetic 
circular dichroism demonstrate a strict coupling of SO and 
CO, present also above Tn and likely even above Tco, 
where diffuse scattering experiments indicate that within 
individual bilayers SO and CO correlations are still strong. 
The presence of coupled SO and CO correlations within in- 
dividual bilayers is a robust feature of the compound. How- 
ever, oxygen off-stoichiometry quickly destroys the inter- 
bilayer spin- and charge-correlations, leading to "glassy" 
behavior 

While charge and spin order have been established, 
they are not really understood: this would imply elucidat- 
ing what drives these orders, whereas currently no theo- 
retical works have considered the experimental charge or 
spin order What is clear is that the charge order cannot 
be driven entirely or even predominantly by electrostatic 
repulsion, because the realized pattern is not the one min- 
imizing this repulsion. Lattice effects as well as magnetic 
exchange appear to be relevant, but the determination of 
the main driving forces remains a challenge principally 
for future theoretical work, though experimental work on 
phonons and magnetic excitations could make an impor- 
tant contribution: this is the most critical issue for future 
investigations. 

Rare earth ferrites other than LuFe20i have been stud- 
ied much less, which is largely due to the absence of single 
crystals that are stoichiometric enough to reveal 3D CO 
and SO. With this hurdle being overcome fTi], progress 
can be expected. CO and SO correlations leading to diffuse 
scattering very similar as in off-stoichiometric LuFe204 
are commonly observed, and for YbFe204 also Mossbauer 
spectroscopy leads to the expectation of overall behav- 
ior very similar to LuFe204. However, for rare earths (or 
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In) with larger deviations of ionic size from Lu, experi- 
ments on high-quality polycrystalline samples indicate that 
completely different CO and SO is established. For exam- 
ple, electron diffraction on YFe204 demonstrates a num- 
ber of competing very complex CO phases that have yet 
to be solved. A partial alternative to single-crystal work 
is to rely on experimental techniques not requiring long- 
range CO and SO, such as XMCD, which has been suc- 
cessfully used to study ion-size effects on spin-charge cou- 
pling 11161 . The intercalated compounds R2Fe^07 have 
been even less studied, but given indications of charge or- 
der stable up to much higher temperatures 111031 deserve 
more attention. A similar tuning as with intercalation may 
be possible with thin films \ 120 121 1 22II123L on which, 
however, 3D charge ordering has yet to be demonstrated. 

Ferroelectricity from charge ordering? While the 
coupled complex charge and spin orders on a geometri- 
cally frustrated lattice are also of high intrinsic interest, 
it is not the property that has attracted most research to 
the rare earth ferrites. The title question "Ferroelectricity 
from iron valence ordering in rare earth ferrites?" accord- 
ing to current experimental knowledge has to be answered 
in the negative, at least for the "prototype" LuFe204. The 
charge order leading to the (i i^)-type reflections does not 
involve polar bilayers, and according to recent dielectric 
spectroscopy studies, the large remanent polarization de- 
duced from pyroelectric current measurements may likely 
be affected by contact effects. Future theoretical work 
leading to an understanding of what drives the experimen- 
tally indicated non-polar charge order would be important 
for a final resolution and to provide hints for other families 
of compounds. In rare earth ferrites other than LuFe204 
the charge order can be different and has not been solved, 
therefore in principle the possibility of ferroelectricity 
from charge ordering remains within the family. How- 
ever, with dielectric behavior quite common among other 
i?Fe204 and even the most likely ferroelectric behavior in 
a doped intercalated compound suggesting a polarization 
three orders of magnitude weaker than what is expected 
for ferroelectricity originating from charge ordering, there 
are no strong indications for success. 

Looking beyond rare earth ferrites, the findings re- 
viewed here provide some general lessons for CO-induced 
ferroelectricity. First and most important, a non-negligible 
residual conductivity is unavoidable in CO materials, par- 
ticularly at and above Tco- Apart from being detrimental 
with regards to applications, this also renders the standard 
macroscopic probes of ferroelectricity doubtful as pointed 
out in t45J and for rare earth ferrites discussed in Sec. 
[3] Therefore it is imperative to obtain also microscopic 
proof of at least a polar crystal structure, better yet of a 
polar crystal structure that can be switched by an elec- 
tric field — though for such in-electric field measurements 
residual conductivity and associated self-heating effects 
are also problematic, see Sec. 14.31 Second, CO materials 
have a propensity for variable stoichiometry, which can 



make the elucidation of the CO crystal structure very dif- 
ficult. Therefore, paying due respect to sample quality is 
crucial. Finally, in CO materials there is often a strong cou- 
pling of charge with orbital and spin degrees of freedom 
leading to complexity that hampers theoretical treatments. 
For rare earth ferrites, theoretical work suggests and/or as- 
sumes ferroelectric CO (see Sec. 14.4) . in contrast to the 
experimental findings. For other materials, there are many 
more examples of compounds proposed to exhibit ferro- 
electricity from CO based on calculations than examples 
of experimental indications of the same. 

With the prototypical example of rare earth ferrites ap- 
parently being a non-example, a valid "proof-of -principle" 
example of a material exhibiting multiferroicity driven by 
charge ordering has yet to be demonstrated experimentally. 
Among the few possibilities in oxides with at least some 
experimental indications most promising seems classical 
magnetite, which does seem to have a charge-ordered crys- 



tal structure that is polar 111241 . although it appears to be of 
a relaxor-type with switching only possible at low T and 
high frequency 01251 . Apart from oxides, low-dimensional 
molecular compounds (see e.g. 11261 ) are increasingly be- 
ing scrutinized for multiferroicity originating from CO. A 
recently proposed 111271 . but also disputed [[1281 . example 
is a metal-organic charge-transfer salt with apparent small 
polarization but no demonstrated polar crystal structure. 
Although the temperature scale in all of the above is far be- 
low room temperature, an unambiguous proof of the mech- 
anism would be important, opening the route to search for 
such materials closer to applicability. 
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